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Abstract 

The classical literature on optimal liquidation, rooted in Almgren-Chriss models, tackles 
the optimal liquidation problem using a trade-off between market impact and price risk. 
It answers the general question of optimal scheduling but the very question of the actual 
way to proceed with liquidation is rarely dealt with. Our model, that incorporates both 
price risk and non-execution risk, is an attempt to tackle this question using limit orders. 
The very general framework we propose to model liquidation with limit orders generalizes 
existing ones in two ways. We consider a risk-averse agent whereas the model of Bayraktar 
and Ludkovski [13j only tackles the case of a risk-neutral one. We consider very general 
functional forms for the execution process intensity, whereas [25] is restricted to exponential 
intensity. Eventually, we link the execution cost function of Almgren-Chriss models to the 
intensity function in our model, providing then a way to see Almgren-Chriss models as a 
limit of ours. 



1 Introduction 

Since the late nineties and the first papers on the impact of execution costs on trading 
strategies {e.g. [T3]), an important literature has developed to tackle the problem of optimal 
liquidation. This literature, often rooted in the seminal papers by Almgren and Chriss [71[H], 
has long been characterized by a trade-off between, on the one hand, market impact that 
encourages to trade slowly and, on the other hand, price risk that provides an incentive to 
trade fast. 

The first family of models, following Almgren and Chriss, considered general instanta- 
neous price impact (sometimes called execution cost) and linear permanent price impact. 
Several generalizations have been proposed such as an extension to random execution costs 
[To] , or stochastic volatility and stochastic liquidity [SI [6]. Also, many objective functions 
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to design optimal strategies have been proposed and discussed in order to understand the 
assumptions under which optimal strategies are deterministic (as opposed to adaptive) . The 
initial mean-variance framework has been expressed in an expected utility framework (us- 
ing a CARA utility function) in [36] and [27], a mean-quadratic- variation framework has 
been considered [1,9^ 48j, initial-time mean- variance criterion has been discussed OUT], and 
the very interesting case of a general utility function has recently been considered in |45| 
to justify aggressive- in-the-money or passive- in-the-money strategies. Slightly different ap- 
proaches have been proposed in this first generation of models (see e.g. ^T] and [33] )• They 
all derive from the initial models by Almgren and Chriss since market impact is either per- 
manent or instantaneous. In other words, they do not take into account the resilience of the 
underlying order book. 

Another family of models appeared following a paper by Obizhaeva and Wang [1^. In 
these models, the limit order book is directly modeled and the authors consider its resilient 
dynamic after each trade. This second generation of optimal liquidation models, based on 
transient market impact, has developed in recent years ([1], [2], [3] and [12]). It raises the 
theoretical question of the functional forms for the transient market impact that are com- 
patible with the absence of price manipulation (see [1], [21] and |23j). 

All these models only make use of market orders, and hence only consider liquidity- 
taking strategies. They do not consider the possible use of limit orders that provide liquid- 
ity, nor the possible use of dark pools. Notwithstanding the preceding criticism, models a 
la Almgren-Chriss provide a rather acceptable answer to the macroscopic question of the 
optimal scheduling of liquidation - at least once the instantaneous market impact function 
has been replaced by an execution cost function modeling the ability to trade over short pe- 
riods of time, with all possible means including limit orders, dark pools and market orders. 
However, they do not answer the question of the optimal way to proceed in practice and 
the methods currently used in the industry are seldom based on optimal control models at 
the microscopic level. This paper provides such a model of optimal liquidation using limit 
orders, and can be used, either to liquidate a portfolio as a whole over a few hours, or on 
shorter periods of time to follow a trading curve, be it a TWAP curve, a VWAP curve or an 
Almgren-Chriss (Implementation Shortfall) trading curve. 

In our approach, a trader post limit orders (thus providing liquidity instead of taking it) 
and does not know when his orders are going to be executed, if at all. As a consequence, the 
classical trade-off between market impact / execution cost and price risk is not central in our 
model. In our setting, a new risk is borne by the trader because execution is now a random 
process. This non-execution risk is very different, in its nature, from price risk. This new risk 
characterizes the recent literature on optimal liquidation, which focuses on the optimal way 
to liquidate rather than on optimal scheduling. The recent literature on optimal liquidation 
focuses indeed on alternatives to the use of market orders. Kratz and Schoneborn |37j pro- 
posed an approach inspired from models of the first family, but with both market orders and 
access to dark pools. Although they considered no risk-aversion with respect to the new risk 
borne by the trader, their model is one of the first in this new family of models. The optimal 
split of large orders across liquidity pools has then been studied by Laruelle, Lehalle and 
Pages in |38j. Liquidation with limit orders has been developed by Bayraktar and Ludkovski 
[13] for general intensity functions but only in a risk-neutral framework. Gueant, Lehalle 
and Fernandez- Tapia |25j considered in parallel the specific case of an exponential intensity 
for a risk-averse agent. More recently, Huitema [M] considered liquidation involving market 
orders and limit orders, and Guilbaud and Pham [30] also proposed a liquidation model in 
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a pro-rata microstructure. 

One should also note that many models dealing with high-frequency market making have 
been developed that can be adapted to deal with optimal liquidation. Building on the model 
proposed by Ho and Stoll [32] and then modified by Avellaneda and Stoikov [ll]0 Cartea, 
Jaimungal and Ricci [18] considered a model with exponential intensity, market impact on 
the limit order book, adverse selection effects and predictable a. Cartea and Jaimungal 
|17] recently used a similar model to introduce risk measures for high-frequency trading. 
Earlier, the same authors proposed a model jl6j in which the reference price is modeled by 
a Hidden Markov Model. Eventually, Guilbaud and Pham [29] also used a model including 
both market orders and limit orders at best (and next to best) bid and ask together with 
stochastic spreads. As it is shown in appendix B, our model can be used to model trading 
on both sides of the market. Our choice to focus on optimal liquidation is mainly justified 
by practitioners' needs. 

In this paper, we generalize both [13] and [25]. We indeed consider both general shapes 
for the intensity functions, and an investor with a CARA utility function. Moreover, we 
present a limiting case in which the size of the orders tends to and we show that this lim- 
iting case is intrinsically linked to the usual continuous framework of Almgren and Chriss, 
although the latter framework only considers market orders. This limiting case helps to 
understand the meaning of intensity functions for quotes corresponding to marketable limit 
orders. 

In Section 2, we present the setting of the model and the main hypotheses on execution. 
The third section is devoted to solving the partial differential equations arising from the 
control problem. Then, in Section 4, we provide illustrations of the model and we exhibit 
the asymptotic behavior of the quotes, generalizing therefore a result presented in [25]. Sec- 
tion 5 is dedicated to the study of a limit regime that corresponds to orders of small size. 
This fifth section leads to results linked to those obtained for the fluid limit in |13] , here in 
a risk-averse setting. This result is exploited in Section 6 that draws parallels between our 
model and the usual Almgren-Chriss framework. 



2 Optimal execution with limit orders: the model 
2.1 Setup of the model 

Let us fix a probability space (fi, T, P) equipped with a filtration {J-'t)t>o satisfying the usual 
conditions. We assume that all random variables and stochastic processes are defined on 

(Ji)t>o,P). 

We consider a trader who has to liquidate a portfolio containing a quantity qo > oi a 
given stock. We suppose that the reference price of the stock (that can be considered the 
first bid quote for example) follows a Brownian motion with a drift: 

dSt = iidt + crdWt. 

To model limit orders and the execution process, we first introduce the set of admissible 
strategies: 

A = {{St)t&[o,T]\{^t)t predictable process ,5" £ L~(ri x [0,r])} . 

See [35] for the solution of the Avellaneda-Stoikov equations. 
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The trader under consideration continuously proposes an ask quote Sf = St + 5t. He 
will sell shares according to the rate of arrival of liquidity-taking orders at the price he quotes. 



His inventory , that is the number of shares he holds, evolves according to the following 
dynamics: 

dqf = -AdN^, 

where A'^'^ is a point process giving the number of executed orders, each order being of size A 
- we suppose that A is a fraction of qo. The intensity process {\t)t of the point process N^, 
that is the arrival rate of liquidity-taking orders, depends on both the (ask) price quoted by 
the trader and the size of its orders: 

At = Aa(S'" - St)lgS_^Q = AA(5t)lg«_>o, 
where A a : M —> IR+ satisfies the following assumptions: 

• Aa is strictly decreasing - the cheaper the order price, the faster it will be executed, 

• lim5_^+oo Aa((5) = 0, 

• Aa G 



As a consequence of his trades, the trader's cash account has the following dynamics: 



dXf = {St + 6t)AdN°. 

Now, coming to the liquidation problem, the trader has a time horizon T to liquidate his 
shares and his goal is to optimize the expected utility of his P&L at time T. We focus on 
CARA utility functions so that the trader considers the following optimization problem: 



supE 



- exp (-7 (X^ + 4(5t - iiq^T)))) 



where 7 > is the absolute risk aversion parameter characterizing the trader, where Xj, is 
the amount of cash at time T and where q^ is the remaining quantity of shares at time T. 
In this setting, the trader can sell the shares remaining at time T in his portfolio at a price 
below the reference price, namely St — (■{qT)^ the function t being a positive and increasing 
penalization function, measuring execution cost. 

We associate to this stochastic control problem the value function Va defined by: 

VA{t, X, q, s) = sup E [- exp f-7 f X^^'^ + - f (4^''^))) ) 

where A{t) is the set of predictable processes on [t,T], bounded from below and where: 

dSl'' = ndT + adWr, S^'' = s, 

dq'f^' = -AdNr, ql^"^' = q, 
the point process having stochastic intensity (At-)t with A,- = AA(6T-)lq^_>o- 

This setting deserves several remarks. First, orders are of constant size A. This is a 
modeling choice corresponding to the way practitioners proceed with liquidation, A being 
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then a fraction of the average trade size (ATS). Also, we imphcitly assume that our orders 
are either entirely filled or not filled at all. In other words, there is no partial fill in this 
model. This hypothesis is a questionable one since partial fills are common in practice. 
When using the model in practice, one can always consider a convex combination of optimal 
quotes between two multiples of A. Allowing for partial fills would make the model more 
realistic. However, it is complicated. 

A second important point regards negative 6s. We indeed assume that Aa is defined on the 
entire real line and not only on M-|-. Our model allows to post orders at a price below the 
reference price. If the reference price is the first bid quote, these orders correspond to mar- 
ketable limit orders. One may then wonder why there is execution uncertainty associated 
to these orders. An answer is linked to high-frequency traders whose capacity to rapidly 
cancel trades forces practitioners to use fill and/or kill orders or other types of marketable 
limit orders and not market orders. Also, considering the entire real line allows to introduce 
indirectly execution costs for liquidity-taking orders. We shall see in Section 6 that there 
is a link between the execution cost functions of Almgren-Chriss models and the intensity 
functions on {5 < 0}. It is noteworthy that if one wants to avoid negative 5, adding a hard 
constraint (5 > is also possible and does not raise any difficulty. This constrained frame- 
work is discussed later in this article (see Section 3.3). To avoid negative 6, some authors 
(see for instance [13] ) considered an intensity function that blows up at 5 = 0. A natural 
consequence of this modeling choice is that there is unlimited liquidity available at 5 = 0. 
This is not a correct approach in our view. 

The third and last point regards the structural assumption Aa(5)A'^(J) < 2A'/\(6)'^. This 
hypothesis, already present in [13], is a sufficient condition to guarantee uniqueness of the 
optimal trading quote. To understand the intuition, let us consider the expected PnL when 
posting an order at a distance 6 from the reference price. This expected PnL is propor- 
tional to 5A/\{6): 6 is the premium over the reference price and A/\{6) is the instantaneous 
probability that a trade takes place at a distance 5 from the reference price. A natural 
condition for this expression to have a unique maximizer is: AaA^ < 2A'a. In our case, the 
inequality can be binding because of risk- aversion. It is noteworthy that our framework can 
be used even if this hypothesis is relaxecH. However, the consequence is that there may be 
multiple optimal quotes. For the sake of exposition, we chose to present the model under 
this structural assumption. 

2.2 A system of ODEs for the value function 

The optimization problem set up in the preceding paragraphs can be solved using classical 
Bellman tools. To this purpose, we introduce the Hamilton-Jacobi-Bellman equation asso- 
ciated to the optimization problem, where the unknown ua is going to be equal to the value 
function Va defined above: 

(HJB) = dtUA{t,x,q,s) + fj.dsUA(t,x,q,s) + ^a'^dggUA{t,x,q, s) 
-|-supAa(s'^ - s) [uAit,x + As",g - A,s) - u/^{t,x,q,s)] , 
with the final condition: 

ua{T, X, q,s) = - exp (-7 {x + q{s - £{q)))) , 
and the boundary condition: 

UA.{t, X, 0, s) = — exp (— 7x) . 
^The same is true for the assumption Aa G C^. 
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Since we use a CARA function, we can factor out the Mark-to-Market (MtM) value x+qs 
of the portfoho. This remark leads to considering the change of variables u/^{t,x,q, s) = 
— exp (— 7(0; + gs + 0A(i, ?)))■ In that case, the above HJB equation with 4 variables is 
(formally) reduced to the following system of ODEs indexed by q: 

(HJe^) = jOtOAit, q) + llJiq - -1 (T q + -ff A ( 1 , 

with 

0A(r,(?) = -%)(?, eA(i,o) = o, 

where 

Ha{p) = supAa((5) (1 - e-^^(^-f)) . 

3 Solving the optimal control problem 

This section aims at solving the optimal control problem set up in the preceding section. 
We first concentrate on the equation (HJ^^). Then, we provide a verification theorem that 
indeed gives a solution to the control problem and characterizes in a simple way the optimal 
quotes. The last subsection is dedicated to the addition of a hard constraint 5 > (J"*'". 

3.1 A solution to (HJgi^) 

We start with a lemma about the hamiltonian function Ha- 
Lemma 1. Let us define La{p,6) = Aa{S) (l - e-T^(^-P)). 

Vp G R, 5 I-)- La{p,S) is strictly increasing on (— oc,5^(p)] and strictly decreasing on 
[S'^{p) , +00) , where ^a(p) uniquely characterized by ^Eg^^ : 

Slip) - 4 log f 1 - 7A^^^^!^) = p, (EsA . 
Moreover, p i-)- is a function. 

Subsequently, Ha is a function with: 

Aa(^1(p))^ 



Ha{p) = 7A- 



Proof: 

Straightforwardly, S ^ La{p,S) is strictly increasing on (— oo,p]. 

Regarding the behavior of this function on [p, +00), let us notice that La{p,p) = and 
that lim5_j._|_oo La{p, S) = 0. 

Now, if we differentiate, we get: 

dsLAip,S) = A'^iS) (1 - e-^^(^-f)) +7AAA(5)e-^^(^-f). 
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Hence, dsL^{p^p) = 7AAa(p) > and there is at least one 5* € {p,+oo) such that 
dsLA{p,S*) = 0. 

Such a S* must satisfy: 



Now, f(x) = X — log ^1 — 7A ^f l^j ^ defines a strictly increasing function 



since 



fix) = 1 + 



A^(x)^-Aa(x)A:^(x) 
^ A'^(x)2-7AAA(x)A^(a;) 

-7AAa(x)A'^(x) JA^(x)2 - Aa(x)A:^(x) 



A;^(x)2-7AAa(x)A^(x) A^(x)2-7AAa(x)A'^(x) 

is strictly positive because of the hypotheses on Aa- 

Hence 6*, defined by f(6*) = p, is unique and La{p, •) is strictly increasing on (— oo, S'^ip)] 
and strictly decreasing on [(5^(p), +oo), where is uniquely characterized by: 

imp)) = hiP) - ^ log fi - 7A^^^SM)) = p. 

Using the implicit function theorem, this also gives that p S'^{p) is a function. 

Plugging the relation for d^ip) definition of Ha then gives the last part of the 

lemma. □ 

Now, we are going to prove a comparison principle for the system of ODEs. This result 
is useful in two ways. First, it gives a priori bounds that will allow us to prove the existence 
of a solution to (HJ^i^). Second, it will provide bounds to ^a independently of A in Section 
5, when we shall consider the limiting behavior of ^a as A ^ 0. 

Proposition 1 (Comparison principle). Let r G [0, T), G N. 

Let 6_A '■ [t, T] X {0, a, ... , fcA} ^ R 6e a function with respect to time with 
VgG{0,A,...,A;A}, dA{T,q)<-l{q)q G [r,r], ^A(i,0)<0, 

and 

< 7Ma(*, q) + im- q +Ha\ 

Let 9a '■ [t, T] X {0, a, ... , kA} be a function with respect to time with 
VgG {0,A,...,fcA}, 9A{T,q)>-i{q)q Vi G [r,T], ^A(i,0)>0, 

and _ _ 

u l222^f. f 9Ait,q)-eAit,q-A) \ 
> jOtOAit, q) + inq - ^1 (T q + a I 1 • 



Then 

0A > 9 A- 

Proof: 
Let a > 0. 

Let us consider a point (tajQa) such, that 

Mtl, q*a) - ^a(C, Qa) - «(r - C) = sup e^{t, q) - 6 Ait, q) - a{T - t). 

{t,q)e[T,T]x{0,...,kA} 

IHIj^T and 7^ then: 

dtedtl,q:)-dteA{tl,q*J<-a. 
By definition of since q^ ^ 0: 

CQ - Mtl, q*a) > iA(.t*a, g: - A) - 9Ait*a,q*a ' A), 

and 

0A{t*a, Q*a) - OA{t*a,q*a " A) > ^a(C, " ^a(C, Qa " A). 

Now, by definition of the functions 9 a and ^a, we have: 

0<l[dteA{tl,q*J-dteA{tl,q*J] 
^ '^^ (^ 9Aitl,q*a)-eA{t*a,q*a-^) ^ _ AK^ql) -6 Aitl^ql - l^) 

Since Ha is a decreasing function, we have 

( lA{tl,ql)-e_A{tl,ql-A) \^ ^ ^A{tl,ql)-eA{tl,ql-A) 

This leads to < —7a which is not possible. 

Therefore = T or g* = 0, so that: 

sup 9A{t,q)-9A{t,q)-a{T-t) 

{t,q}e[T,T]x{0,...,kA} 

= max( sup 9 a{T, q) -9 AiT, q), sup 9 A{t,0) -9A{t,0) - a{T -t)] <0. 
\qe{Q,...,kA} telr,T] J 

Thus V(t, q),9^{t, q) — 9A{t, q) < a{T — t)< aT. Sending a to proves our result. □ 

We are now ready to prove that the equation (HJ^^) has a unique solution. 

Proposition 2. There exists a unique function 9a ■ [0, T] x AN — )■ M such that: 

• t E [0,T] t-^ (^a("j ?))geAN is continuously differentiahle. 

• 9a is a solution of (HJ^t^ ) . 
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Proof: 



We proceed by induction on q. For q = 0, we have by definition 9A{t, 0) = 0. 

Now, for a given q G AN*, let us suppose that 9a{-,(i') '■ [0,T] ^ M is a function 
yq' <q- A. Then, the ODE 

n u ^^ ^ 2 22^0- f eA{t,q) - 0A{t,q - A) 

= jdt9A{t, q) + ^nq - T;Tcr q +Ha' 



2 ' ^ " V A 
with the terminal condition 

9A{T,q) = -£{q)q, 

satisfies the assumptions of Cauchy-Lipschitz theorem. Consequently, there exists a unique 
solution t 9A{t,q) on a maximal interval that is a sub- interval of [0,T] and we want to 
show that this sub-interval is [0, T] itself. 

To prove this, let suppose by contradiction that {i, T] is the maximal interval with i >0. 
Let us notice that, because Ha is positive, t i->- 9A{t,q) — lJ.q{T — t) + ^ja'^q^{T — t) is 
decreasing. Hence, the only possibility for {i, T] to be a maximal interval in [0, T] is that 
limj_^j+ 9 Ait, q) = +00. 

Now let us consider > 0, A; = ^ and r = i +r). We define on [r, T] x {A, ... ,q} the 
two functions 9^ and 9a defined by: 

9a = 9a 

and 

W < q,9A{t, q') = IJ,+q{T - t) + -//A(0)(r - t). 

7 

These two functions satisfy the assumptions of the above comparison principle. We 
indeed have that: 

Vg' < q,9AiT,q')=0> -eiq')q' G [r, T], ^a(<, 0) = fi+q{T - 1) + -Ha{0)(T - 1) > 0, 

7 

and 

, r.-^ , / 1 2 2 /2 , fOA{t,q')-9Ait,q'-A)\ 
Vt G [T,T],yq G {A, . . . ,q},'ydt9Ait,q)+-flJ,q - -1 cf q +Ha I 1 

= -7-^a(0) - 7M+9 + ll^q' - \l^(y^q^ + Ha{^) 
7 ^ 

= -lii^'^q - m') - ^7^0-^?'^ < 0. 

Hence, Vry > 0, Vi G [i + r/, T],9A{t, q) < IJi^q{T + I^Ha{^){T - i), in contradiction with 
the fact that \mi^^i+ 9 Ait, q) = +00. 

Hence, t 0a(^, q) is defined on [0, T] and this proves the result. □ 
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3.2 Verification theorem and optimal quotes 

Now, we can solve the initial optimal control problem and find the optimal quotes at which 
the trader should post his limit orders. 

Theorem 1 (Verification theorem and optimal quotes). Let us consider the solution Oj\ of 
the system (HJg^). 

Let us define UA{t, x, q,s) = — exp(— 7(x + qs + 9/\{t, q))). 



We have: 

• UA is a solution to (HJB), 

• UA is equal to the value function Va ■ 

Moreover, the optimal ask quote = St + S'^{t), for qt > 0, is characterized by: 



Slit) = SI 



f OA{t,qt)-eAit,qt-A) 



V A 

where is the function defined in Lemma 1. 

Proof: 



From the very definition of and ua, it is straightforward to see that ua is a solution 
of (HJB). 

We indeed have that the boundary condition and the terminal condition are satisfied, 
and for g > A, we have that: 

dtUA{t,x,q,s) + iJ,dsUA{t,x,q,s) + ^a'^dggUAit,x,q, s) 

+ supAa((5) [uAit,x + As + AS,q — A,s) - UA{t,x,q,s)] 
S 



1 



= -7UA{t, x, q, s)dt9A{t, q) - jfiquAit, x, q,s) + -a q UA{t, x, q, s) 

+ sup AA{S)uA{t, X, q, s) (exp {-^{Ad - {9 Ait, q) - 9 Ait, q - A)))) - 1) 
5 



= -UAit,x,q,s) 



ldt9Ait, q) + - 2 ^^7^9^ 



+ sup Aa(<5) (1 - exp (-7(A(5 - (0A(t,g) - ^A(t,g - A))))) 
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-UAit,x,q,s) 



u \^ 1 2 2 2 ^^. (9Ait,q)-9Ait,q- A) 

ldt9Ait, q)+im- 1 Q. +Ha' 



0. 



Now, we need to verify that ua is indeed the value function associated to the problem 
and to prove that our candidate ((5^)t is indeed the optimal control. To that purpose, let us 
consider a control S G Ait) and let us consider the following processes for r G [t,T]: 



dSi'' = iidT + adWr, S^'" = s. 



t.s 



^j^t,x,S ^ ^ dr)AdNr, X, 



t.x.S 



X, 
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t,q,5 

Qt =1^ 



where the point process has stochastic intensity (A,-),- with A,- = AA((5r)lq^_>oIl 
Now, let us write Ito's formula for uaEI 



ua{T, x!^^^, q'rj?^ , S^rf') = UA{t, X, q, s) 



+ jj""" (nA(r, X*'^'^ + AS'/ + - A, S'/) - u^ir, X'/J , q';^/ , S'/)] 



where M is the compensated process associated to N for the intensity process (At-)t- 

Now, we have to ensure that the last two integrals consist of martingales so that their 
mean is 0. To that purpose, let us notice that dsU = —'jqu, and hence, since the process 
q^''^'^ takes values between and q, we just have to prove that: 



E 



E 



and 



1^ UA(.T,X'/^',q'/J,S'/fdT <+oo, 
j^^ \uAir, Xfj^ + AS'f + A5,, q'f_l^ -A,S: 

£\uA{T,X';^^',q'/J,St'n 



Xrdr 



< +00, 



X-rdr 



< +00. 



We have: 

ua{t, X'/'',q'/'', S'/f < exp (27||0a||oo) exp (-2-f{X'/'' + q'/'^S'/) 



< exp (27||0a||cx>) exp i^-2jix - q\\6 lU + 2*7 ^mf^^ S'/l.^^^^^^^^^ ^t,.^^; 

< exp (27||0a||oo) exp (-27(x - q\\S-\\oo)) (^1 + exp (^-2^q ^mf^^ 5*'^^^ 



Hence: 



E 



ruA{r,Xt'-'',qi!i'\Sin'dT] =E [ ^ UAir.X'/J.q'^'^J ^S'ffdr 
Jt J VJt 



<exp(27||0A||oo)exp(-27(x-^||r||oo))(r-t) 1 + E 



exp -27g inf S ' 



t,s 



< +00, 



because of the law of inf^g[( -pj St'^ . 
■^Tliis intensity being bounded since 6 is bounded from below. 

""^The equality is still valid when Qt = because of the boundary condition for ua, and because the intensity 



process is then assumed to be 0. 
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Now, the same argument works for the second and third integrals, noticing that 6 is 
bounded from below and that A is bounded. 



Hence, since we have, by constructiorH 



we obtain that 



E 



/rp •yt,X,5 t,q,5 rit,, 



: E 



UA{T,Xr^J ,qrft ,6^ ) < UA{t,x,q,s), 

and this is true for all 5 € A{t). Since for 6t = S^{t) we have an equality in the above 
inequality by construction of the function 5^, we obtain that: 



sup E 



ua{T,x!^^'^ ,qll?'^ ^St/f) <UA{t,x,q,s)=W. ua{T, X!^'^'"'^ ,q'rf^'"'^ , S^r^'^) 



I.e. 



sup E 



exp(-7(X^^'' + g^'^''(5^^-%^^'')) 



< UA{t, X, q,s)=E [- exp (-7 (^X^f"'^'^ + q^^'^'^iS'/ - l{q%^'^*^)) 
This proves that ua is the value function and that t 1— )• 6*^{t) is optimal. 



□ 



Theorem 1 proves that the optimal quotes are deterministic. This is linked to the use of 
a CARA utility function as in the usual Almgren-Chriss framework. Theorem 1 also pro- 
vides a simple way to compute the optimal quotes. One has indeed to solve the triangular 
system of ODEs (HJg^) to obtain the function 9a- Numerically, this does not constitute 
any difficulty and one may use for instance a Euler scheme. Then, once has been com- 
puted, the optimal quotes are given by the simple expression 5\ where 

the function 5*^ is implicitly characterized by the equation i^s*^ of Lemma 1, and can be 
easily computed using Newton's method for instance. 



3.3 Introducing a hard constraint 5 > (5™*" 

In the above framework, 6 was allowed to take any value on the real line. To avoid mar- 
ketable limit orders, one might want to impose a constraint 6 > (J*"*" where J*"*" would be 
positive. Using the same tools as above, the problem with the additional constraint 5 > S"^^"" 
can be solved easily. 



The (HJB) equation becomes 

= re, g, s) + x, q, s) + ^^^^^^^-(t, x, q, s) 



^This inequality is also true when the portfolio is empty because of the boundary conditions. 
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+ sup Aa(s"-s) [^^^"(^,^ + As^g-A,s)-nS*'^(^,a;,g,s)] , 

with the final condition: 

<-(T, X, q,s) = - exp (-7 {x + q{s - iiq)))) , 
and the boundary condition: 

u'^'^it, X, 0, s) = - exp (-7x) . 

We consider the change of variables u^^"'{t, x, q,s) = — exp (—7(0; + qs + q))) , as 

above. We then obtain the following system of ODEs indexed by q: 



(HJ™) = -fdtOr'it, q) + -f^iq - ^-^^a\^ + H^'' 



01'''{t,q) - ef-'^{t,q - A) 



with 



where 



et''{T,q) = -l{q)q, ^r"(i,0)=0, 
i/r"(p)= sup AA((5)(l-e-^^(''-^')) = sup La(p,(5). 



The important point here is to recall that 5 ^ La{p, S) is strictly increasing on (—00, S\{p)] 
and strictly decreasing on +00). Hence, the unique maximizer of (5 i-> L^{p,5) over 

{5 > J'"^'^} is max(5'"*'^,5^(p)). 

Let us define p"^^" = (5"*^" - log (l - 7A^jJ|^) . The hamiltonian function iJ^"'" 
can be written: 



^ ^ ' |Aa(<5™") (1 - e-^^('5"''"-f)) ifp<p^ 



mm 



mm 



It is a locally Lipschitz and decreasing function. In particular, the counterpart of Propo- 
sition 1 and Proposition 2 holds: there exists a unique function 0^"" solution of (HJ^mm). 
Therefore, we can enounce a verification theorem and find the optimal quotes. The proof is 
mutatis mutandis the same as for Theorem 1. 

Theorem 2 (Verification theorem and optimal quotes). Let us consider the solution 0^*" 

of the system (HJ 

Then: 

- exp(-7(a; + qs + O^'it, q))) = sup E [- exp (-7 fe"''' + (^j?^\s'^' - ^(g^^'^))))! , 

where A^'^^(t) is the set of predictable processes on \t,T], hounded from below by (J"*'" and 
where: 

dSl'' = ndT + adWr, 5^*'" = s, 
^Xt,x,s ^ ^ Sr)AdNr, Xl'""'^ = x, 

the point process N having stochastic intensity {Xt)t with A,- = A./^{Sr)'i-q^_>o- 
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Moreover, the optimal ask quote Sf = St + (5^*"*(t), for qt > 0, is characterized by: 
5T^*{t) = max ( 5^'^,~5l " " 



where ) is the function defined in Lemma 1. 

4 Examples and properties 

4.1 The case of an exponential intensity function 

In the above section, we generalized a model already used in [25J, in which the intensity 
functions had exponential shape: Aa((J) = A^e~^^^ . 

In the case of exponential intensity, we can write the results of [25J (in a slightly more 
general case than in the original paper) in the language of this paper. In fact, the reason why 
closed-form solutions can be obtained in the exponential case is that the equation (HJg^) 

simplifies to a linear system of equations when we replace the unknown 6/^ by exp ( 



Proposition 3. Assume that h.^{5) = A/^e 



Then, Ha{p) = ^{i + ^) ' ''^^A^e-^^P. 

Also, if we consider 9^ : [0, T] x AN — )• M, the unique solution of (HJg^), then 
= exp ^^^A^ is the unique solution of: 

dtWA{t,q) = --^kAfJ'qwA{t,q) + ^7^Acr^g^w^A(t, g) - ^A ( 1 + j WA{t,q- A), 

with 

^«A(T,g) = e-^^('?)^ WA{t,Q) = l, 
and the optimal quote, for qt > 0, is given by: 

X* m 1 1 f^!!A(t^qt)\ ^ 1 . A ^ 7A\ 
KA \WA[t,qt-^)J 7^ V «^A/ 

4.2 Numerical examples 

We now provide numerical approximations of both the function 9 Ait, q) and the optimal con- 
trol function 5*^{t,q). These numerical approximations allow to compare what happens in 
the pure exponential case and what happens when another intensity function is considered, 
especially for negative 5s. In this section, we consider as an alternative to the exponential 
form for Aa a functional form A. a (see Figure [I]) that prevents the use of marketable limit 
orders. The intensity function Aa prevents the use of marketable limit orders since Ka is 
constant (in fact, decreasing very slowly to satisfy the hypotheses of the paper) for negative 
5. Also, we included the commonly observed fact that the probability to be executed does 
not correspond to the exponential intensity framework for small positive 5s. 
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Intensity function 

0-14 -I 




0.02 H — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 
-1 1 2 3 4 5 

delta (in lick) 



Figure 1: Intensity functions (here represented on [—1,5]). Line: Aa(5) = Ae A = 0.1 (s 
k = 0.3 (Tick-^). Dotted line: AaI*^), identical to Aa(5) for S>2. 




50 100 150 200 250 300 50 100 150 200 250 300 



Figure 2: Solution 0A(t,g) for Aa (top) and Aa (down), go = 400, A = 50, T = 300 (s), ^ = 0, 
a = 0.3 (Tick.s-^), 7 = 0.001 (Tick"^) and i{q) = i = 3 (Tick). The index q E {0, 50, . . . ,400} 
of each curve can be read from the terminal values. 



Figure [2] and Figure [3] represent respectively the solution 6a and the optimal quotes 
^A^t^Q) as given by Theorem 1^ From Figure [H we know that ^a is not a monotonic 
function of q. It is important here to recall the economic meaning of 6a- The certainty 
equivalent of holding q shares at time t is qs + 6A{t,q). Hence, 6A(t,q) is a risk-adjusted 
value of holding q shares at time t in excess of the MtM value qs. The reason why 6A{t,q) 
is not a monotonic function of q can then be understood easily. At the time horizon T, the 
function is decreasing but far from T two effects are at stake. On the one hand, when there 
are many shares in the portfolio, there will be many trades and hence more opportunities 
to make money through limit orders: this goes in the direction of an increasing function 
6 Ait, •)• On the other hand, the larger the inventory to liquidate, the more price risk. This 
goes in the direction of a decreasing function 6A{t, •) since it is a risk-adjusted value. 



^One may wonder why we choose a risk aversion parameter 7 = 
fact an important impact since the shares are sold by groups of 50. 
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0.001. This figure seems small but it has in 




50 100 150 ZOO 250 300 50 100 150 200 250 300 

Time (sec) Time (sec) 



Figure 3: Solution 6l{t,q) for Aa (top) and Aa (down), qq = 400, A = 50, T = 300 (s), = 0, 
cr = 0.3 (Tick.s~2), 7 = 0.001 (Tick~^) and £{q) = i = 3 (Tick). The lower the quotes, the higher 
q E {50,..., 400}. 

Although, there is almost no difference between the two cases Aa and Aa as far as ^a is 
concerned, this is not true anymore when it comes to the optimal quotes 5a- We indeed see 
on Figure [3l as expected, that, in the case of the intensity function Aa, there is no negative 
optimal quotes. This very conservative choice for the intensity function is a way to avoid 
any influence of the intensity function on the set {5 < 0}. Also, we see that, since Aa is not 
of exponential form for small positive 6s, the lower bound for the optimal quotes is higher 
than expected. 

Another way to prevent marketable limit orders is to impose 6 > (J™*"- = as in Section 
2.3. In that case, if we consider the intensity function Aa and the same parameters as above, 
we obtain the quotes given on Figure HI These quotes are almost exactly the same as if we 
had floored the optimal quote 6X of the unconstrained problem to 0. 



Opiimal quotes 




50 100 150 200 250 300 



Time (sec) 

Figure 4: Optimal quote for Aa when the constraint 5 > is imposed, qo = 400, A = 50, 
T = 300 (s), ^ = 0, a = 0.3 (Tick.s"^), 7 = 0.001 (Tick"^) and £{q) = i = 3 (Tick). The lower 
the quotes, the higher q G {50, . . . , 400}. 
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4.3 Asymptotic quote 

In [25], we obtained a limiting regime when T — )• cx). This result generalizes to our general 
framework. More exactly, we obtain: 

Proposition 4 (Asymptotic behavior). Let us suppose that: 

• limp^+oo-f^A(p) = 0. 

• 7>0, 

• ^ < ^-/a^A. 

Then, the asymptotic behavior of is: 

hm 9^(0, q)=A Yl ^a' f ^T^^g'^ - jf,q') = 9'S{q). 

g'e{A,2A,...,g} ^ ^ 

The resulting asymptotic behavior of the optimal quote is: 

^lirn^<5X(t = 0)=6l (^H^^ {^I'^'ll " 7^go)) = S*^ ■ 

Proof- 
Let us define for q £ AN: 

g'G{A,2A,...,g} ^ ^ 

Let us define for t > and q € AN, 9^, the unique solutiorj^] of: 

n P,rn 1 2 2 2 ^^ f 9lit, q) - 9l{t, q - A) 

= -7dt9A{t,q) +jfiq- --y a q + a ' 



2 ' ^ ^ V A 
with 

Ol{0,q) = -£{q)q, 0X(t,O)=O. 
Then, because we just reversed time, we want to prove that: 

VgeAN, lim 9Ut,q) = e^{q). 

We proceed by induction. The result is true for q = 0. Let us suppose that the result is 
true for q — A for some q G AN* . 
Then: 

Ve > 0, 3ig_A, Vt > tg_A, |0a(*> 9 - A) - e^{q - A)| < e. 
Since H/\ is a strictly decreasing function, we obtain that Vt > tg-A^ 

1 2 2 2 ^ f 9Ut,q)-9^{q-A) + e 



''This is guaranteed if lim^^+oo (5Aa(<5) — 0. 

^This is well-defined because of the assumptions on Ha and on the parameters, and because 7 > 0. In the 
risk-neutral case, there is no upper bound for the optimal quotes when T — > 00. This constitutes an important 
difference between our model and the model of Bayraktar and Ludkovski |I3j . 

^To prove that this function is well-defined, one can use the same tools as in Proposition [21 
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or equivalently: 





6^{q)-9^{q-A) ^ 
^ f0^iq)-e^{q-A) 



<ldtei{t,q)<HA 




Hence, Vt > tq^^'- 



e'^{t,q)>e^{q) + e^dte^^{t,q)<0, 



and 



9'^{t,q)<e^{q)-e^dte'^{t,q)>0. 



As a consequence, if there exists t' > tq^^ such that \9'^^{t',q) — 0'^{q)\ < e then. 



^t>t',\e\(t,q)-Q^{q)\<e. 

In particular, if \Q\(tq_A,q) - Q'^{q)\ < e then, Vt > tg__A, \0'^^{t,q) - 0'^{q)\ < e. 



Now, if 6^^{tq-A,q) > (9) ~^ ^' tl^sn there are two possibihties. The first one is that 
the function t > tq-A '-^ ^a(*i9) decreasing and in that case it is bounded from below by 
0^{q) — e and must converge. Since limt_>.+oo ^a(*' q~^) — ^a (^~^)' ^he only possible limit 
for 9^^{t,q) is 6'^{q). The second possibility is that t > tq-A ^ ^a(^'9) * decreasing 

function and in that case there must exists t' > tq-A such that ^a(^'' i) — "I" ^- Since 

6'a(*''9) > ^Ail) - we now obtain that Vt > tq, |6l^(t,g) - 9'^{q)\ < e. 

Finally, if 6^(tq-A, q) < — e, then there are two possibilities. The first one is that 

the function t > tg_A '-^ ^a(^' i) increasing and in that case it is bounded from above by 
9 Ail) + f S'^d must converge. Since limt^+00 9^j^{t,q — A) = — A), the only possible 
limit for 6'^^(t,q) is 9^{q). The second possibility is that t > tq^A ^ 9''^{t,q) is not an in- 
creasing function and in that case there must exists t' > tq^A such that 9^^{t', q) > 9'^{q) — e. 
Since 6IX(t', q) < 9'^{q) + e, we now obtain that Vt > tq, |6l^(t, q) - e'^{q)\ < e. 

The conclusion is that lim supj^_,_o<3 \6^{t,q) — 9'^{q)\ < e. Sending e to 0, we get the 
result for ^a- 

The result for the optimal quote is then straightforward. □ 

These asymptotic formulae deserve some comments. Firstly, regarding the above discus- 
sion on monotonicity, we know that Ha is a decreasing function, therefore, the asymptotic 
limit ^A (■) either a decreasing function (when ^7^0"^ A^ — 7/iA > Ha{0)) or a function 
that is first increasing and then decreasing (otherwise) - see Figure [5] in our case. Secondly, 
coming to the optimal quotes and the role of the parameters, we can analyze the way 
depends on /i, a, 7 and Aa. The best way to proceed is to use the expression for Ha found 
in Lemma 1 and to notice that an equivalent way to define 5'^ is through the following 
implicit characterization: 



t^qo + ^icr'^qo = A 



Aa{S*^? 



7AAA(<5r)-A'z 
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Figure 5: 0^{q) for Aa defined above, go = 800, A = 50, yU = 0, a = 0.3 (Tick.s 2), 7 = 
0.001 (Tick"^). 

It is then straightforward to see that 5^° is an increasing function of fi. A trader ex- 
pecting the stock price to go up is indeed encouraged to slow down the liquidation process. 
Similarly, we see that 5'^ decreases as a increases. An increase in a corresponds to an in- 
crease in price risk and this provides the trader with an incentive to speed up the execution 
process. Therefore, it is natural that the asymptotic quote be a decreasing function of a. 
Differentiating the above expression with respect to 7, we see that the asymptotic quote 
decreases as the risk aversion increases. An increase in risk aversion forces indeed the trader 
to reduce both non-execution risk and price risk and this leads to posting orders with lower 
prices. 

Now, if one replaces the intensity function A a by AAa where A > 1, then it results in an 
increase in 5'^ . This is natural because when the rate of arrival of liquidity-taking orders 
increases, the trader is more likely to liquidate his shares faster and posting deeper into the 
book allows for larger profits. 



4.4 The influence of A 

In addition to the asymptotic regime, we can consider different sizes A of orders. 

We see on Figure [6l that there is little difference between the two cases we considered. 
This is linked to the existence of a limit regime as A — > and the next section is dedicated to 
its analysis. The limiting equation for will turn out to be a classical equation in optimal 
liquidation theory (see Section 5). 



5 Limit regime A — > 

In the preceding sections, the size of the order posted by the trader was constant equal to 
A, a size that is supposed to be small with respect to qq. As a consequence, the question of 
the limiting behavior when A tends to is relevantly To that purpose we need to make an 
assumption on the behavior of the intensity function with respect to the order size A. The 
"right" scaling (already used above for the numerics underlying Figure [6]) is to suppose that 



'Although this is not recalled, it is assumed that A is always chosen as a fraction of qg. 
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Figure 6: Optimal quotes for q E {50, 100, . . . ,400}, for go = 600, T = 1200 (s), /i = 0, a = 
0.3 (Tick.s-^), A = 0.1 (s-i), k = 0.3 (Tick-^), 7 = 0.001 (Tick"^) and i{q) = i = 3 (Tick). 
Line: A{6) = Ae'^^ and A = 50. Dotted line: A{5) = 2^6"'=^ and A = 25. 



With this scahng , along with additional technical hypotheses, our goal is to prove the 
following Theorem that is rather technical and echoes the results obtained by , here with 
risk aversion (i.e. 7 > 0) whereas [13] deals with the risk- neutral case: 

Theorem 3 (Limit regime A ^ 0). Let us suppose that: 

• A{5)A"{5) < 2A'{5f 

• lim^^+oo 5A{6) = 

• i is a continuous function 

For a given A > 0, let us define 9'^ on [0, T] x [0, go] f^V- 



that is 



'0A(t,O), ifq = 0, 

fAit,ik + l)A), ifq(£{kA,{k + l)A]. 

Then 6'^ converges uniformly toward a continuous function 9 : [0, T] x [0, go] - 
the unique viscosity solution of the equation (HJum): 

-jdt9{t, g) - 7M + h^a^q^ - H {dg9{t, g)) = 0, on [0, T) x (0, go], 
9{t,q) = 0, on[0,T]x{0}, , (HJu^) 

9{t,q) = -l{q)q, on {T} x [0, go], 

where H{p) = 7SUP5 A(5){6—p) and where the terminal condition and the boundary condition 
are in fact satisfied is the classical sense. 

To prove this theorem, we first need to study H and the convergence of the Hamiltonian 
functions towards H. We start with a counterpart of Lemma 1 that requires A{6)A"{6) < 
2A'{5f. 
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Lemma 2. Let us define L{p,d) = A{6) {6 — p). 

Vp G M, (5 H-)- L{p,5) attains its maximum at S*{p) uniquely characterized by: 

^ ' A'{5*ip)) 
Moreover, p i-^ 6* {p) is a function. 

Subsequently, H is a function with: 

-A'(5*(p)) 

Proof: The proof is similar to the proof of Lemma 1. 

Now, we can state a result about convergence that also provides a uniform bound for the 
hamiltonian functions: 

Lemma 3. H/^ converges locally uniformly towards H when A — >■ with G R, i?A(p) < 
H{p). 

Proof: 

1 —At 

For a fixed x > 0, the function / : A G i-> — is a decreasing function (/(O) = x). 
Hence, 

1 _ p-7A(<5-p) 1 _ -7A'(<5-p) 

VO < A' < A, sup K{5) < sup K{6) < 7 sup A(5) {6-p). 

5>p ^ S>p ^ S>p 

This gives: 

Ha{p)<Ha'{p)<H{p). 

Now, because H is continuous, using Dini's theorem, if we prove that convergence is 
pointwise, convergence will be locally uniform. We then only need to prove pointwise con- 
vergence of i?A toward H. Using Lemma 1 and Lemma 2, we see that it is sufficient to prove 
that (5^ converges pointwise towards 6*. 

For that purpose, notice that the sequence of functions f/^{x) = x — ^^ log ^1 — 7A-^^^ 
is an increasing sequence of increasing functions. Hence, by the unique characterizations of 
{p) ^* {p) ' ^^^^ {P) increases as A decreases to and is bounded from above 

by 5*{p) which is the only possible limit. Hence 5*^{p) ^ as A — >■ and this proves 

the result. □ 

Now, we provide a uniform bound for the 9^ that will be important in the proof of 
Theorem 3. 

Proposition 5 (Bounds for 6'a). G [0,T],Vg G {0, A, . . . ,^0}, 

-^(go)9o - /^~go(r -t)- ^7f^'9o(r -t)< e^it, q) < fi+qo{T - t) + -H{0){T - t). 

2 7 
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Proof: 

To prove these inequalities, we use the comparison principle of Proposition [TJ 
If^A(t,(?) =/i+go(T-t) + ii/A(0)(r-t), then: 

^a(T,(?) = 0> OAit,0)=fj.+qo{T-t) + -H^{0){T-t)>0, 
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1 2 2 2 ^rr ( /^{t, q) - 9 ^{t, q - ^) 

idtOAit, q) + -ffiq - -7 0"^ + Ha ' 



and 

^r)tf) /\ (t. n\ A- ^ii.n — , „ „ , , , . 

2 ' ^ V ^ 

= -7-^a(0) - Jfi^qo + im - \l^(y^q^ + ^^a(0) 
7 2 

= -7(m^9o - /^g) - ^7^cr^9^ < 0. 
Therefore, ^aC*, q) < 0A{t, q) = /U+<7o(T - t) + i/fA(0)(r - t). 
The uniform upper bound is then obtained using Lemma 3. 
Now, if ^A = -£{qo)qo-l^-qo{T-t) - ^-fa^q^{T - t), them 
O-AiT, q) = -Kqo)qo < -i{q)q, O^it, 0) = -£{qo)qo - f^~qo{T -t)- ]p^a^ql{T - t) < 0, 



f^a \ , 1 2 2 2 , TJ ( lA{t,q) - l^it^q- ^) 

ldtO^[t, q) +lliq- -1 CF q +Ha' 



and 

jdtO^{t,q)+-fnq - ^, - . . , ^ 

- |-'-222| l222it7 fn\ 

= 7M qo + ^1 qo + 7M9 - 2^ q + Ha{0) 

> j{fi-qo + ^iq) + \l^(y^{ql - q^) > 0. 
Therefore, 9A{t, q) > ^Ait, q) = -£{qo)qo - f^'qoiT - t) - \l(J^ql{T -t). □ 
We are now ready to start the proof of Theorem 3. 
Proof of Theorem 3: 

We first introduce the fohowing half-relaxed limit functions: 

9{t, q) = lim sup lim sup sup | 0a(*'i I*' ~ *l + I?' ~ Q'l ^ ~ f > 

e(t,q) = lim inf lim inf inf 1(9^ \t' - t\ + \q' - q\ < - \ ■ 
j-s>+oo A-i>0 [ j } 

6 and are respectively upper semi-continuous and lower semi-continuous and the goal of 
the proof is to show that they are equal to one another and solution of the partial differential 
equation (HJiim)- 
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Step 1: 9 and 9_ are respectively viscosity subsolution and viscosity supersolution of the 
equation (HJiim)- 

To prove this point, let us consider (i*, (?*) and a test function (p such that 6 — (p attains a 
local maximum at {t* , q*). Without loss of generality, we can assume that (p{t* ,q*) = 0{t* , q*) 
and consider r > such that: 

• the maximum is global on the ball of radius r centered in (t* ,q*). 

• outside of this ball, (p > 2sup^ ||^a||oo ~ this value being finite because of the uniform 
bound obtained in the above Proposition. 



Following Barles-Souganidis methodology [12], we know that there exists a sequence 
{An,tn,qn)n such that: 

• A„ ^ 0, {tn,qn) {t*,q*) 

• 0''^^{tn,qn)^0{t*,q*) 

• 6'^^ — (j) has a global maximum at {tn,qn). 

Now, because of the definition of 9'^ , if {t*,q*) G [0,T) x (0, go] we can always suppose 
that qn > A„ and tn ^ T and then, by definition of 6a„- 

-lOteAAtn, qn) - imn + Ip qn - ^An ( 1=0. 

Because -ffA„ is decreasing we have, by definition of {tn,qn)' 



■ ^A„fa,gn) - 6'A„fa,gn - An) \ ^ f ct){tn, qn) - ^tn, qn - ^n)^ 



Similarly, since tn < T , we have, for h sufficiently small: 

Oa„ {tn + h, qn) - 6a„ (tn, qn) < 4>{tn + h, qn) - 0(t nj qn)- 

Hence: 

dt9A,A'tn,qn) < dt(j){tn,qn)- 

These inequalities give: 

o , /. X 1 2 2 2 u [ (t){tn,qn) - (t>itn,qn- An)\ 
--idt(p{tn,qn) - imn + ^ q^ - Ha„ ( 1 < 0. 

Using now the convergence of {tn,qn) towards {t*,q*) and the local uniform convergence 
of Ha.„ towards H, we eventually obtain the desired inequality: 

-jdtcP{t*,q*) - + ll^a^q*^ - H {d,cP{t\q*)) < 0. 

We see that the boundaries corresponding to q = qo and t = play no role. However, we 
need to consider the cases t* = T and q* = 0. 

If t* = T and q* then there are two cases. If there are infinitely many indices n 
such that tn <T then the preceding proof still works. Otherwise, for all n sufficiently large, 
GAn{tn,qn) = —^{qn)qn and hence, passing to the limit, 6{t*,q*) = —i{q*)q*- 
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Eventually, we indeed have that: 



mm{-jdtct>{T, q*) - 7/xg* + -j^a\*^ - H {d^4>{T, q*)) , e{T, q*) + £{q*)q*) < 0. 

If g* = then there are also two cases. If there are infinitely many indices n such 
that > and t„ < T then the initial proof still works. Otherwise, for n sufficiently 
large 0^^^ {tn,qn) = —£{qn)Qn or OA„{tmQn) = and hence, passing to the limit, we obtain 
e{t*,q*)=0. 

Eventually, we indeed have that: 

mm{--fdtcP{t*,0) - H {dq4>{t\Q)) ,9it*,0)) < 0. 

We have proved that ^? is a subsolution of the equation in the viscosity and one can 
similarly prove that ^ is a supersolution. 

Step 2: Vg G [0, qo] , ^(T, q) = e{T, q) = -i{q)q. 

We consider the test function (/){t,q) = Ce{T — t) + ^{q — qref)^ where g^e/ ^ [0, 5o] is 
fixed, where e > is a constant and where Q is a constant that depends on e, and that will 
be fixed later. 

Let {te, qe) be a maximum point of ^ — ^ on [0, T] x [0, go]- Then: 

e{T, qref) < e{te, qe) - C,{T - Q ' ^{Qc ' Qref? ■ 

This inequality gives q^ qref as e — >^ 0. 
Now, 

-ldt(t){te,qe)-llJ.qe + ]p'^(^'^qe^-H{dq(t){te,qe)) = ^Ce-^ IJ-qe + ^p"^ (t"^ Qe^ - H {^{Qe - Qref)^ 

>jCe-H(-^-^)+ M (-J^^q+lJ'aV). 

Hence, if Ce is positive and greater than (^—^^ — inf^gp^g,,] (—7^9 + f 7^o"^g^) + 1, we 

see that we must have either te = T and 9{te, qe) < —^{Qe)Qe or (for sufficiently small e, only 
in the case where qref = 0), 9{te,qe) < 0. 

If qref 7^ 0, we then write, for e sufficiently small: 

e{T,qref) < 0{t„qe) - Ce{T - te) - ^{qe " Qref)^ < 0{te,qe) < -i{qe)qe. 

Sending e to we obtain 

0{T,qref) < -iiqref)qref- 

If qref = 0, then we know from the above two inequalities that: 

e(T, qref) < eite, qe) " CeiT - te) " ^(^6 - Qref)^ < 0{te, qe) < max(-£(ge)(?e, 0) = 0. 
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Hence, Vg G [0, qQ\^6[T, q) < —£{q)q and the same proof works for the supersolutior0 to 
get yq€[0,qo],e{T,q) > -£{q)q . 

As a consequence Vg € [0,qo],6{T,q) < —i{q)q < 0{T,q), and eventuahy, because 
0{T,q) < 6{T,q), we obtain that G [0, go],^(r, g) = -i{q)q = 9{T,q). 

Step 3: 'it e [0,T], e{t,0) = 

Concerning the boundary condition corresponding to g = we can apply the same ideas 
but only to the supersolution 9. 

Let us consider indeed t^e/ £ [0)^) and the test function <j>{t,q) = —C^q — \{t — trefY • 
Then, let (t^, q^) be a minimum point of — (/> on [0, T] x [0, q^]. We have: 

This inequality gives t^ — )• t^ef as e — ?> 0. 
Now, 

-^dt4>{te, qe) - 7We + Ii^^^q' - H {dg(t>{te, Qe)) = 2j^-^^ - ^fiq, + ^7^2^,2 _ H i-a) 

< 27^ - 7;U(?o + ^7''^''?o' - i-Ce) . 

Since liuip^^ao H{p) = +oo, we can always choose > so that the above expression is 
strictly negative. 

As a consequence, for e sufficiently small, we must have = and ^(t^, 0) > 0. 
Consequently: 

0{tref, 0) > e{t„ q,) + C,q, + \t, - tref? > 0. 

This result being already true for tref = T, we have that 0{t,0) > 0,Vt G [0,T] and in 
fact 9{t, 0) = 0, Vt G [0, T] because of the definition of e{t, 0). 

Step 4: Vt G [0, T], there exists a sequence {tn,qn)n such that t„ 7^ i, 7^ 0, {tn,qn) 
(t,0), and 0t„,g„) ^0 

To prove this claim, we prove that QA^t) = O^iT — t, A) converges uniformly (in t) toward 
as A 0. 

By definition, 5a(0) = -£(A)A and g'^{t) = //A - i7a2A2 + ^H^ (^) • 
We now distinguish two cases: 
Case 1: /x < 0. 



^The only difference is that we have to pass to the limit in the case qref = to obtain the conclusion. 
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The stationary state of the above ODE is = Ai7^"'^(^7^(T^A^ — 7/iA) and g/s, is in- 
creasing on {^A < S-A }• Since, i:^H~^^ {^-^"^ a"^ /^^ - 7/iA) > l:^H^ {l^f"^ a"^ l:^^ - 7//A) as soon 
as A < A', g"^ is positive for A sufficiently small. As a consequence, since 5a(0) < 0, g^ is 
increasing on [0,T]. 

Now, g'^{t) < lH(^^y Hence, 

and this gives g/^{t) < AG"^ (^^^ , where G{x) = j^dy. 
Consequently, 

VtG[0,T], -i{A)A<g^it)<AG-'(^-^y 

Since lim^^+oo 6A{6) = 0, wc have linia^^+oo H(x) = 0, and therefore liuix^+oo = +00. 
This proves that gA converges uniformly toward on [0, T]. 
Case 2: /it > 0. 

In this second case, we know that BAq such that A € (0, Aq) t— > /iA — ^^/a^A'^ is an 
increasing and positive function. Now, let us introduce Ai G (0, Aq). VA € (0, Ai), we have: 



dy 



7/xAi- 17^2^2 + 
9A(r) 



-dy = T. 



-eiA)A 7/^A-i7VW + /fA(i) 

This gives 5A(r) < AG^^ (:^), where Gai(x) = ^^Ai-i7^a^ Af 
Now, because lim3;_^+oo H{x) = 0, we have: 

liminf^^> ' 



Hence, 



a;->+oo x 7/xAi — i72cr2A^ 



T 1 

liminf ^ > 

A-)-0 



This gives 
and sending Ai to 0, we get: 



7AGk: (^)-7/^A,-i7V^Af 



limsup5A(r) < r(//Ai - Iw^Aj), 



limsup5(A(r) < 0. 
A-^-o 



The result is then proved since: 



vte[o,r], -^(A)A<(7A(i) <5A(r). 

Step 5: Comparison principle: 6 < 6. 
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Let us consider a > and the maximum M = max(( q)g[o,r]x[o,go] ^(^' Q)~^i^i q) — a{T—t). 
If m = maxjg[o,T] ^{'t, 0) — 6{t, 0) — a{T — t) < M then we distinguish two cases. 

Case 1: m < 0. 

We introduce ^,{t, q, t' , q') = ^(t, q) - e{t' , q') - a{T - t) - - ^^=|^. 



Let us consider {tc,qe,t'^,qi) a maximum point of We have M < ^^{te,qe,t[,q'^) and 
we are going to prove that hm'mi^^Q ^^{te, qe,t'^, q'^) < 0. 

We have ^^{te,q^,t^,q^) < ^eite,Qe,te,qe) and hence iSizpl + ih^pl jg bounded. As a 
consequence, — t'^ ^ and q^ — q'^ ^ 0. 

Now, if for aU e sufficiently smaU we have (i^, g^, t^, G [0, T) x (0, go] x [0)^) x (Oj^o] 
then we have: 

— ,-'-2 2 2 zrr rile ~ 9e \ ^ o 

c ^ c 

and 

-27 £ - 7/xgg + ^6 - ^(2 -) > 0. 

Hence 7a — 7/x(ge — q'^) + \'y^a'^{q^ — q'^) < and this is a contradiction as we send e to 

0. 

The consequence is that there exists a sequence ^ such that (ig^, ge„) *1„) I'en) verifies 
tg„ = r,Vn or = 0, Vra or 4„'Vra or gg^,Vn. 

Then 

Km sup $g„ (tg„ , gg„ , i'g„ , J < Km sup ^(tg„ , ^g J - ^(t'g„ , J - a(r - tg„ ) 

n— >'+oo n— ^+00 

< max e{t, q) - OU, q) - a(T - t) < maxfO, m) < 0. 

Hence in that case M < 0. 
Case 2: m > 0. 

In that case, we replace ^ by ^ + m in the above case, and we obtain, instead of M < 0, 
the inequality M <m which contradicts our hypothesis. 

It remains to consider the case m = M. We know then that the maximum M is attained 
at a point {tmax,^) and we suppose that M > 0. 

Prom Step 4, we consider a sequence {tn,qn) such that 7^ tmax^ Qn 0) (tn^Qn) 
{tmax,0) and e{tn, qn) 0. 

We define: 

^nit, q, t', q') = e{t, q) - 0, q') - a{T - t) - j^^P' I " ' ^ 
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This function attains its maximum at a point (i* , g*, t'* , g'* ). We first consider tlie 
inequality 'i>n{tmax,0,tn,qmax) < ^n{tn,Qn^t'*„,Q'n)- 

d{tmax, 0) — 0{tn, Qn) — Ol{T — tmax) — \ in tffiax I 

(f* _ /'*12 / /* _ * 

I ^max \ \ Qn 

We then have t* - t'* ^ and q'n - Qn ^ ^■ 

Hence limsup„ , g*) — 9{t''^,Q'n) ~ ~ O ^ Now, the maximum M is also 
given by lim„ 9{tmax,^) — i.{tn, Qn) — (^(T — tmax), and we obtain that the penalization term 

^" ^" - 1 converge to 0. 



qn 

This gives g'* = g * + gn + o(g„) > 0. 

Now, if we have infinitely many n such that t'* = T, then: 

M = lime{tmax, 0) - e{tn, Qn) - "(T - tmax) < SUp e{T, q) - 9{T, q) = 0, 
" <?G[0,go] 

and this contradicts our hypothesis. 
Otherwise, for all n sufficiently large: 

-^^jf'-rn - ^^^n + h'^'r: -h(-- -A)>^- 

l^-n ^max] ^ \ Qn \ Qn J J 

Now, going to the subsolution, if there are infinitely many n such that t* = T, then: 

M = \mi'd{tmax,^)-9itn,Qn)-Oi{T -tmax)< SUp e{T, q) - 9{T, q) = 0, 
" qe[0,qo] 

in contradiction with our hypothesis. 

Else, if g* = and 9{t*^.q*^) < for infinitely many n, then we obtain M = lim„ 9{tmax,0) — 
9{tn,Qn)~C({T—tmax) — \tn—tmax\ < straightforwardly, and this contradicts our hypothesis. 
Hence, the viscosity inequality must be satisfied and we get: 

I In ''max | \ 9n \ Qn J J 

Combining the two inequalities eventually leads to 07 < as n — )■ oo and this is a con- 
tradiction. 

We have obtained that M < and hence 9 — 9< aT. Sending a to 0, we get 9 < 9. 

This proves that ^ = ^ is in fact a continuous function that we call 9, solution of the 
PDE (HJiiin). Using the same techniques as above, it is clear that 9 is the unique viscosity 
solution of (HJ iim ). 



We have: 



9(t,q) = 9(t,q) < liminf 6'^(i,g) < limsup6'^(i, g) < 9(t,q) = 9(t,q). 
A->0 A^O 
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Hence 0{t,q) = lim^^Q 9'^{t,q) and, by the same token, liniA-s>o,(t',g')^(i,g) ^a(*'' ~ 
0{t,q) so that the convergence is locally uniform and then uniform on the compact set 
[0,T]x[0,qo]. □ 



6 Link with Almgren-Chriss 

6.1 Interpretation of the PDE in the Umit regime 

In the above section we proved that 9^ converged to 9, which is the unique continuous vis- 
cosity solution of (HJiim), the limit condition and boundary condition being satisfied in the 
classical sense. 

Now, we are going to link this equation to a classical equation of Almgren-Chriss-like 
models. The intuition behind the link between our framework in the limit regime A — > 
and the Almgren-Chriss framework is that non-execution risk vanishes as A tends to 0. 
Hence, the only remaining risk is price risk, corresponding to the term ^j'^a'^q'^ in the above 
equation. To see more precisely the correspondence between the two approaches, let us write 
the hamiltonian function as: 



H{p) = 7supA(5)(5 -p) 

s 

= 7supv(A~^(i;) — p) 
= supvA~^{v) — pv, 

where the last equality holds since lim5_!.+oo SA{5) = 0. 

Hence, if we define for v > 0, f{v) = — A~^(t;), then we can define the function H{p) = 
sup^>Q —f{v)v — pv and write the partial differential equation for 9 as: 

' -dt9{t, q)-tiq+ IjaV " H {dq9{t, q)) = 0, on [0, T) x (0, q^], 
9{t,q) = 0, on[0,r]x{0}, 
[9it,q) = -£{q)q, on {T} x [0, go]- 

This equation is the Hamilton-Jacobi equation associated to the Almgren-Chriss optimal 
liquidation problem with an instantaneous market impact function (per share) / and with 
a final discount i{qT) per share. More precisely, the above Hamilton-Jacobi equation is the 
Hamilton-Jacobi equation associated to the optimization problem^ 



inf r (-q'(t)f{-q'{t)) - f,qit) + lw\itf) dt + q{TY{q{T)), 



geAC{0,T),q{0)=qo Jq 

where AC{0,T) is the set of absolutely continuous functions on [0, T]. 



However, the instantaneous market impact function / has here a rather unusual form 
since f{v) is negative for v < A(0) and positive for v > A(0), whereas it is usually a positive 



^^In the usual Almgren-Chriss framework, liquidation is mandatory but the theory can easily be adapted to 
allow for a penalization term. 
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function. This must be interpreted in a very simple way: if one needs to obtain an instan- 
taneous volume lesser than A(0), then one will choose a positive 6. In other words, he will 
post a classical limit order, since we assumed that the reference price is the first bid quote 
- this makes sense since non-execution risk disappears in the limit regime A — )• 0. On the 
contrary, if one needs an instantaneous volume greater than A(0), then one will rely on a 
marketable limit order (<5 < 0). 

The above discussion only makes sense at the limit, when non-execution risk does not 
exist anymore. However, it clarifies the meaning of negative 6s. In particular, the above 
correspondence between our model and a model d la Almgren-Chriss provides a possible way 
to solve one of the main practical problems of the model discussed in [25| : the interpretation 
of the intensity functions for negative values of 5. 

6.2 Discussion on the choice of A 

The model we discuss in this paper does not consider explicitly market orders or limit orders 
but rather considers that there is, for each price s° = s + (5, an instantaneous probability to 
obtain a trade at that price. In practice, this interpretation is perfectly suited to classical 
limit orders, but we need to provide an interpretation for the intensity function A on the 
entire real line. In practice, since the model has been designed to liquidate a position with 
limit orders, it should not be used if the liquidation evidently requires liquidity-taking orders 
from the very beginning. However, it may happeUjbecause of a slow execution, that the 
optimal quote in the model turns out to be negativqlfl after some time. 

This issue of negative 6 was present in the model with exponential intensity functions 
introduced in [25]. Although the exponential form was justified for many stocks for positive 
6s, the intensity function was also of exponential form for 6 < 0, and this choice was dictated 
by mathematical needs rather than by empirical rationale. Since A (or in practice Aa) can 
be chosen in our setting, we can improve the initial model. 

First, using statistics on execution, we can estimate the probability to be executed at any 
positive distance from the first bid limit. In practice the profile of the empirical intensity 
for positive 6 is decreasing and may not be convex, especially when the bid-ask spread is 
large (this is the rationale underlying our choice - for 5 > - on Figure [TJ. Then, once the 
function A a has been calibrated for positive 6, several natural choices are possible for Aa((J) 
when 6 < 0. Instead of extending the function for negative 6 using a specific functional 
form as in [25], we can assign to Aa((5) a constant value when 5 < as above for Aa- This 
corresponds to a very conservative choice that basically prevents the use of marketable limit 
orders since, intuitively, the optimal quote will then always be positive. Another choice 
consists in using the parallel made between the usual literature and our framework in the 
limit regime A ^ 0. If we indeed omit non-execution risk, we can consider, for 6 < 0, that 
Aa((5) = ^f~^{—6) = ^ sup{u > 0, f{v) < —6} where v f{v) is an instantaneous market 
impact function (average execution cost per share) that is typically equal to nought for small 
values of v and increasing after a certain threshold. 

This choice for Aa is however subject to several comments. First, the function Aa must 
satisfy the hypotheses of the model. In particular, it must be decreasing. However, it may 
happen that, although the specifications for positive 6 and negative 6 are both decreasing, 
the function is not decreasing on the entire real line. Since the function 6 < i-^ 

Although quotes evolve continuously between execution times, using the model in practice requires to post 
orders and to keep them in the order book for some time. Hence, the optimal quote at some point may in practice 
be strictly negative. 
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can be considered a lower bound to Aa(5) because we have to take into account the risk 
of non-execution, we can always scale the function 5 < i— )• Aa((^) so that the resulting 
function Aa is decreasing (and strictly decreasing if we smooth the function). In general, 
the inequality A"A < 2A'^ may not be satisfied but the model can still be used although the 
optimal quotes 6* may not be unique. Second, a question remains regarding the interpreta- 
tion of the model when an optimal quote 5* turns out to be negative. A possible answer, in 
line with the parallel made with Almgren-Chriss-like models, is to send a market order of 
size A{6*) = AAa((5*) (or in practice to use marketable limit orders to obtain this size). 



Conclusion 

In this paper, we analyze optimal liquidation using limit orders. The classical literature 
on optimal liquidation, following Almgren-Chriss, only considers optimal scheduling and a 
new strand of research has recently emerged that uses either dark pools or limit orders 
to tackle the issue of the actual optimal way to liquidate. Our paper provides a general 
model for optimal liquidation with limit orders and extends both |13] that only considers 
a risk-neutral framework and [25] that was restricted to exponential intensity functions. 
Our general framework also sheds new light on the important topic of negative quotes. An 
important improvement of our model would consist in linking the Brownian motion which 
drives the price and the point process modeling execution. Research in this direction has 
recently been made by Cartea, Jaimungal and Ricci [18] to model market making and it 
may adapt to our case. 

Appendix A: The multi-asset case 

This appendix is devoted to the generalization of our results to the case of a portfolio with 
multiple stocks. Our main result (Theorem 1) generalizes to the multi-asset case. 

We consider a trader who has to liquidate a portfolio made of d different stocks with a 
quantity Qq of stock i (i G {1, . . . , d}). 

We suppose that the reference price of stock i evolves as: 

dSi = n'dt + a'dWi, 

with Y(W^, W^) = ip''^)i<i,j<d definite positive. 

The trader under consideration continuously proposes an ask quote for each stock (5^'" = 
SI + 51 for stock i), and will hence sell shares according to the rate of arrival of liquidity- 
taking orders at the prices he quotes. 

The state of the portfolio is (^q^'^^ , . . . ,q'^'^''^ . It evolves according to the following 
dynamics: 

yi,dqf =-A'Ni'^\ 

where N^'^' is a point process giving the number of executed orders for stock i, each order 
on stock i being of size A* - we suppose as above that A* is a fraction of Qq. The intensity 
process (AJ)t of the point process A^*'^' is given by: 
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where the function A^j : M — > R_|_ satisfies the same assumptions as in the single-stock case. 
We suppose that the point processes N^'^^ , . . . , A'''^''''* are independent. 
The cash account X^^'-'^"^ of the trader has then the following dynamics: 

d 



1=1 



Eventually, the optimization problem is: 

/ / d 



sup E 

(5^,...,5d)(iA-^ 



exp -7 4 + Y.^f{Si, - l\cif)) 



i=l 



where the functions £^,. . . ,£'^ satisfy the same assumptions as in the single-asset case. 

This multi-asset setting deserves two remarks. First, the execution processes associated 
to different stocks are independent. Hence, the only difference between the multi-asset case 
and the single-asset cases has to do with price risk: in this multi-asset framework, optimal 
quotes will depend on the correlation structure between stocks. Second, the trader can only 
sell shares, although buying shares might sometimes reduce price risk in practice. 

The counterpart of Proposition 2 and Theorem 1 in this multi-asset framework is the 
theorem below: 

Theorem A.l (Verification theorem and optimal quotes). There exists a unique solution 
^Ai,...,A<*; in time, of the system: 

VtG[0,r),(g\...,g'^)^(0,...,0), 

d 

= 75teAi,...,A'^(i,g',---,9')+7E/^V-27' E P''WgY 

1=1 l<i,3<d 
, A , ( dA\..Mi, g\ . . . , g^) - e^r_^a{t, q\...,q'-A\...,q'') \ 



i=l 

and the conditions: 



Oa\..Mt^ <i\---^<f) = -Y. ^(«')«'' ^A^.M"^, 0, . . . , 0) = 0, 

1=1 

where 

yi,Hi,{p) = supAX,((5') (l - e-^^'('^'-f)) . 
// /^d is this function, then 



Ua^ A'iit,x,q^,...,q'^,s^,...,s'^) = - exp ( -7 ( x 



is the value function of the optimal control problem, and the optimal ask quotes are charac- 
terized by: 

Vz, [d^ijt = o^i 
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where S^^i{p) is uniquely characterized by: 



7A' 



■log I-7A' 



Proof: 



The proof is mutatis mutandis the same as in the single-stock case. Existence of a local 
solution t !->■ 0j\i Ad{t,q^, . . . ,q'^) comes from Cauchy-Lipschitz. To obtain existence on 

[0,T], first notice that e^i_^4t,q\ . . ■ ,q'')+Eti f^'9KT-t)-hEi<ij<dP'''(^'<^'l'lHT~ 
t) is a decreasing function of t. Now, using a comparison principle similar to Proposition [U 
we have: 

d -] 

e^i_^4t, q\...,q'')<J2 -t) + -Y^ H^,.{0){T - t), 

i=l ^ 1=1 

and this bound guarantees that there is no blow up. This lead to global existence. Unique- 
ness follows from a comparison principle. 

As far as verification is concerned, the proof is exactly the same as in the single-asset case. □ 



Appendix B: Trading on both sides 

Along with optimal liquidation, an important strand of research on optimal trading is high- 
frequency market making. We claimed in the introduction that many models can be used to 
deal with both optimal liquidation and market making. We illustrate this claim and show 
that our framework can easily be adapted to two-sided trading. We present the model in the 
single-stock case. The multi-stock case works the same but notations make the exposition 
cumbersome. 



We consider a stock with a reference price following a Brownian motion with a drift: 

dSt = jidt + adWt- 

The trader under consideration continuously proposes a bid quote = St — 5\ and an 
ask quote S'f = St + 51 . His inventory q^^'^" evolves according to the rate of arrival of 
liquidity-taking orders at the prices he quotes: 

dqf^'^ =^dNf -MNf, 

where iV^' and N^'^ are the point processes giving the number of executed orders respectively 
on the bid side and on the ask side, each order being of size A (on both sides). The intensity 
processes (A^)f and (Af)j of the point processes N^^ and A'"^" are given by: 

A = ^A{St - St)l Sb^s'^ = AA((5t)l 5b,S<^ 



and 



At = ^A{St - St)l = Aa(5")1 Sb,sa 

It- ^ ^ It- ^ 



where Aa : M — )• satisfies the same assumptions as in the case of the optimal liquidation 
model, and where Q is a bound on the inventory. The bounds on the inventory have two 
roles: (i) they stand for the risk limits the traders have in practice, and (ii) they allow to 
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write a verification theorem as in 

As a consequence of his trades, the cash account X^"'^" of the trader has the following 
dynamics: 

dXf = {St + df)AdNf - {St - 5\)MNf. 
Finally, the optimization problem of the high-frequency market maker is: 



sup E 



exp^-7^X^' +qj,' ST-\qT Klkr 1)^ 
where I satisfies the same assumptions as in the case of the optimal liquidation model. 



The counterpart of Proposition 2 and Theorem 1 in this market model framework is the 
theorem below: 

Theorem B.l (Verification theorem and optimal quotes). There exists a unique solution 
Oa, in time, of the system: 

G [0,r),VgG {-Q,...,-A,0,A,...,Q}, Q = ^dt9A{t,q) + J^lq - ^-f^a\^ 

rr f eA{t,q)-eA{t,q + A) \ , ^ f eA{t, q) - eA{t, q - A) 

+ 1<?<Q-HA I I + lq>-QHA I 

and the terminal condition: 

9A{T,q) = -i{\q\)\q\, 
where Ha{p) = sup^ Aa(5) (l - e'^^^^'P'') . 

If 9a is this function, then 

UA{t, X, q,s) = - exp (-7 {x + qs + 9A{t, q))) 

is the value function of the optimal control problem, and the optimal bid and ask quotes are 
characterized by: 



b,, 7* f9A{t,qt)-9A{t,qt + A) 



iS^A)t = S 



A 



A 

9A{t,qt) -9A{t,qt- ^) 



A 

where S'^{p) is uniquely characterized by the equation [Eg*^ ) of Lemma 1. 
Proof: 



The proof is close to the proof in the one-sided case. Existence of a local solution 
t i-> {9A{t, q))q^^_Q _A A Q} comes from Cauchy-Lipschitz. To obtain existence on 
[0, T], first notice that 9A{t,q) + fJ'q{T — t) — ^ja^q'^{T — t) is a decreasing function of t. 
Then, using a comparison principle, we have: 

9A{t,q) < \l2\Q{T-t) + -HA{0){T-t), 

7 

and this bound guarantees that there is no blow up, hence global existence. Uniqueness 
follows from a comparison principle. 

As far as verification is concerned, the proof similar to the proof of Theorem 1. □ 
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